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The spinless resonant level model is studied when it is coupled by hopping to one of the arbitrary 
number of conduction electron channels. The Coulomb interaction acts between the electron on the 
impurity and in the different channels. In case of repulsive or attractive interaction the conduction 
electrons are pushed away or attracted to ease or hinder the hopping by creating unoccupied or 
occupied states, respectively. In the screening of the hopping orthogonality catastrophe plays an 
important role. In equilibrium the weak and strong coupling limit the renormalizations are treated 
by perturbative, numerical and Anderson- Yuval Coulomb gas methods. In case of two leads the 
current due to applied voltage is treated in the weak coupling limit. The presented detailed study 
should help to test other methods suggested for non-equilibrium transport. 
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I. INTRODUCTION 

In the recent years the quantum impurity problem out- 
of-equilibrium has attracted great interest. The most rel- 
evant realizations are the quantum dots connected to at 
least two metallic leads^ and short metallic wires con- 
taining magnetic impurities 2 -. In the impurity problem 
the exact methods play distinguished roles especially the 
Bethe Ansatz and conformal invariance. The generaliza- 
tion of these methods to out-of-equilibrium situations are 
the most challenging new directions. Mehta and Andrei 
are aiming to solve the Kondo problem on a dot with two 
contact attached. First a simple resonant level without 
spin was studied to test the new generalization of the 
Bethe Ansatz method^. Their elegant suggestion is very 
provocative. In order to test this kind of new methods 
we perform a detailed study of that problem using differ- 
ent weak coupling perturbative methods combined with 
NRG. As the final goal we calculate the current flow- 
ing through the impurity when a finite voltage is applied 
on the contacts. The most challenging claim of Mehta 
and Andrei is that the current is a non-monotonic func- 
tion of the strength of the Coulomb coupling between the 
electron on the dot and conduction electrons in the two 
leads. 

In order to make the comparison more explicite we 
generalize the time-ordered scattering formalism for non- 
equilibrium in the next leading logarithmic order. In 
this way the current is calculated as a function of the 
applied voltage and the Coulomb coupling strength. In- 
creasing the Coulomb coupling strength we find also a 
non-monotonic feature but the order of increasing and 
decreasing regions is the opposite to the finding of Mehta 
and Andrei^. 

The model to be treated is the following: A single im- 
purity orbital is coupled to two reservoirs of Fermi gas 
via hopping but the two reservoirs have different chemi- 
cal potentials hl and hr on left and right of the impurity 



in a one dimensional model, lit, — lir — eV is determined 
by the applied voltage V (e is the electronic charge) . The 
Coulomb interaction acts between the electron on the im- 
purity level and the conduction electrons at the impurity 
position. Thus the Hamiltonian has the form 



with 



H = H Q + H 1 +H 2 , (1) 



H = ^2 (k - k a )v F a\ a a ka +e d d^d , (2) 



where k > and k^ — kR = eV/vp, vp is the Fermi veloc- 
ity, a ka is the creation operator of the spinless Fermion 
in lead a — L/R, while Ed is the energy of the local level 
and S is the creation operator for the electron on that 
site. The interaction term is 



H 1 = U(<ftd-±) 



(3) 



>a=L,R 



where U is the Coulomb coupling which in a physical case 
U > 0, a a = -j— ~^2 k aka, and L is the length of the chain. 
The existence of the substraction of 1/2 is not essential, 
they can be omitted and than Ed is shifted as Ed — U/2 
and a local potential — |t7 is acting on the electrons, but 
the latter one can be taken into account by changing the 
electron density of states in the leads at the position of 
the impurity. 

The hybridization between the lead electrons and the 
localized electron is described by 



H 2 = V a ^2 (Jc 



(4) 



where V a is the hybridization matrix element. 

In case of equilibrium it is useful to generalize the 
model to N reservoirs instead of L, R, and then a runs 



2 



through a = 0, 1, . . . , N — 1 and \i a = \i. Then the hy- 
bridization term in i? 2 is chosen in a specific form 

H 2 = V (d^a + aldj , (5) 

indicating that only the electrons with a — are hy- 
bridizing while the others are taking part only in the 
Coulomb screening. Namely, only those electrons are hy- 
bridizing which have the symmetry of the localized or- 
bital (s-like). As a result of the screening the electron 
gas is polarized depending on the occupation of the lo- 
calized state and that polarizations lead to orthogonality 
catastrophe^. 

The model with N = 1 is known as a resonant level 
model and has been studied in great detail^ and the one 
with N > 1 has been introduced to study finite range 
interaction in 3D 7 . 

The goal of the present paper is to provide weak cou- 
pling results for V ^ 0. But before doing that the V = 
equilibrium case is studied in the weak coupling limit 
by diagram technique and then to extend the results 
for stronger couplings Wilson's numerical renormaliza- 
tion group (NRG)£ and Anderson- Yuval Coulomb gas 
method 9 in order to check the validity of weak coupling 
results concerning a specific behavior. Namely at some 
stage of the calculation in the exponent of the renormal- 
ized quantities a combination 

- QoU + \N{e U) 2 (6) 

appears. For U > that is changing sign at QqU = j? 
and that leads in changing the increasing to decreasing 
behavior but that crossover is well beyond the validity of 
the perturbation theory at least for N = 2. 

In order to judge the real situation, an NRG study 
will be performed including the weak (goU -C 1) as well 
as strong coupling regions (qqU > jjj) to get an insight 
whether the crossover indicated above is expected or it 
is likely artifact of the weak coupling theory. We also 
map the problem to a one-dimensional Coulomb model 
closely following the work of Anderson and Yuval, where 
the screening can even be in the strong coupling limit. 
All of these methods are suggesting a coherent picture 
of the crossover and they agree very well especially for 
N = 4. 

The study of such a crossover is especially relevant as 
in the work of Mehta and Andrei^ such a crossover is 
suggested in the current flowing in the non-equilibrium 
case V 7^ at qqU ~ 2. If we could find the crossover 
already in equilibrium then it is obvious to expect that 
in the non-equilibrium situation. 

The paper is organized in the following way: In Sec- 
tion [TT] we provide the analytical perturbative method up 
to next to leading logarithmic order introducing extra 
channels for screening, where the non-monotonic com- 
petion of the vertex and self-energy correction is already 



demonstrated in equilibrium. In Section IIIII the equilib- 
rium calculation is extended to strong coupling by us- 
ing Wilson's numerical renormalization group technique 
and the result is compared to that of the analytical cal- 
culation. In Section IIVI the Anderson- Yuval method is 
presented. In Section [V] the time dependent scattering 
method is applied for non-equilibrium closely following 
the generalized version of Anderson's poor man's scaling 
in the next to leading order and the current is calculated. 
In Section I VII the results are summarized. In Appendix 
A some cancellation due to Ward indetities are discussed. 



II. PERTURBATION THEORY: WEAK 
COUPLING LIMIT 

The resonant level model is given by Eqs. |T]), @ and 
That does not contain non-commuting terms, thus 
Kondo behavior is not expected in the weak coupling 
limit. In the strong coupling limit that model, however, 
can be mapped to an anisotropic Kondo model^i 7 - but 
that mapping is not considered here. The model shows 
strong similarities to the X-ray absorptio n 10 ! 11 , as the 
strength of the interaction (invariant charge) between 
conduction electrons and the electron on the impurity 
level is scale invariant. The system shows scaling in terms 
of reduction of the conduction electron band width, D. 
In case N — 1 the scaling equations were derived by 
Schlottmann 6 and those can be easily extended for arbi- 
trary N. 

There are two different kinds of vertex corrections de- 
picted in FigQ]a,b,c, where the solid lines stand for con- 
duction electrons, the dotted line for electron on the im- 
purity level and the interactions are indicated by dots 
(U) and crosses (V). In case of conduction electrons the 
channel indices are also indicated. The Hartree-Fock en- 
ergy shift can be incorporated by Sd- The self-energy of 
the electron on the impurity is depicted in FigfT] d. In 
the calculation of the self-energy counter terms are intro- 
duced to eliminate the constant terms to keep Sd — un- 
renormalized. Closely following the earlier work o 6 > 12 ' 13 , 
the invariant charge for the Coulomb interaction takes 
the form 

U inv = T{io/D)d(iu/D) , (7) 

where T is the vertex function and 

d{w/D) = G{w 1 D)(u-e d ) 

can be determined perturbatively starting with 1 where 
G is the renormalized one-electron Green's function. The 
functions T(w) and d(u>) are 

r(w) = i + JVegi/ 2 iog(-) + ... 

d(w) = l-Ng 2 Q U 2 log(-) + ... , (8) 

where go is the conduction electron density of states, for 
spinless electrons in one of the channels n = 1, . . . , N — 1. 
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FIG. 1: Vertex diagrams (a,b,c) and the impurity self energy 
(d). Solid lines stand for conduction electron propagators 
while the dotted line for those of electron on the impurity 
level. The interactions are indicated by dots (U) and crosses 
(V). In case of conduction electrons the channel indices are 
also indicated. 



n=0 




n=0 

FIG. 2: Vertex correction to the hybridization. 



As the Coulomb interaction is independent of n, thus 
the factor N occurs. As the consequence of the Ward 
identity relating the vertex correction and the self-energy 
depicted in Figfl}: and d cancel out in {7i nv ^ i 10 ' 1:L i 13 

(T(w)d(w)) = /rt^vH - , (9) 



thus U lnv — U. The renormalization group gives 




for arbitrary n. 

The hybridization contains vertex correction for n = 
(Fig© ' 



V(cu) = V [1 + Ug \og- + ... , (11) 



but it does not contain linear contribution in log(Z)/w) 
in the order of U 2 . The relevant invariant charge Vi nv (u)) 
is 

V inv {to) = V(w)d 1/a (w) (12) 
as the interaction is connected by only one impurity line. 
Thus the terms linear in log(D/u>) are 



1 ; m U 1(1- QoUlog — - In(q U) 2 log - 



(13) 



The result of the renormalization equation is 
V^{u) = v{-) 



(14) 



The second term in the exponent appears as reduction 
of Vi nv (cu) describing the Coulomb screening in the N 
channels. 

For U < the Vmv (w) interaction is always decreasing 
as u) is reduced but for U > it depends on the strength 
of U. For large enough U the screening dominates thus 



decreasing for Uqq < 
Vi nv (lu) = ^ increasing for < Uqq 
decreasing for Uqq > jj 



(15) 



That behavior will be further discussed in Sec lVIl 

The scaling regions, however, are not unlimited as the 

impurity level has its own width, rj mp . There are two 

contributions to the level width rj mp . 

The hybridization broadens the impurity level just in 

case of the Anderson model and that is in the second 

order in V 



r im pM = irg a V 2 (uj) 



-"(b) 



-2 eo U+N(g U) 2 



(16) 

where also the effect of renormalization is taken into ac- 
count. There is also a Korringa like broadening due to 
the creation of electron- hole pairs thus rKoninga ~ U 2 lu, 
where u comes from the phase space. That is impor- 
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FIG. 3: (Color online) The renormalized hybridization as a 
function of qqU for different channel numbers. The interme- 
diate maximum can diverge only for go U — 1, in all the other 
cases the increases are also rather moderate. 



tant only for large to where everything is smooth thus 
the broadening is not effective. 

The broadening due to the hybridization cuts off the 
renormalization procedure at energy u ~ Fj mp (w). Com- 
bining Eqs. (|14p and (|16p and inserting the condition 
given above provides the final V Yon value as 



^ren 

~v~ 



D 



(17) 



As shown in Fig [31 for U < V ICU /V < 1 renormal- 
izes downwards and for even more negative U down to 
zero. For U > and QqU < 2/N first V lcn increases with 
increasing U but for q$U > 2/N it starts to decrease 
and tends to zero again. The intermediate maximum 
appears at QqU = 1/N. For N — 2 that maximum is, 
however, already outside the weak coupling limit, where 
the calculation cannot be trusted. The question still re- 
mains unanswered whether the non-monotonic behavior 
for U > can be traced in strong coupling calculation or 
not. The conclusions for the crossover might be trusted 
only for large N ^> 1, which does not have physical rele- 
vance. 



III. NRG APPROACH FOR V = 

In order to determine the region of validity of the weak 
coupling approach in equilibrium, we have performed nu- 
merical renormalization group£ (NRG) analysis for the 
N = 2 and N = 4 cases. 

In Wilson's NRG technique -after the logarithmic dis- 
cretization of the conduction band- one maps the original 
Hamiltonian of an impurity problem to a semi-infinite 
chain with the impurity at the end of the chain. As a 
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FIG. 4: (Color online) Upper panels: impurity density of 
states for V/D = 0.015 and U/D = 0, . . . , 0.3 as obtained by 
perturbative RG and Wilson's NRG. The lower panel shows 
the renormalized value of the hybridization, Ken as a function 
of the interaction strength U. The numerical data supports 
the weak coupling results for U/D < 0.2. 



consequence of the logarithmic discretization the hop- 
ping amplitude along the chain decreases exponentially 
as t n ~ A - ™/ 2 where A > 1 is a discretization parameter 
(we have used A = 2 throughout the calculations) while 
n is the site index. The separation of energy scales pro- 
vided by the exponentially vanishing hopping amplitude 
allows us to diagonalize the Hamiltonian iteratively to 
approximate the ground state and the excitation spec- 
trum of the full chain. Since we know the eigenenergies 
and eigenvectors of the Hamiltonian, we can calculate 
dynamical quantities such as the density of states using 
the Lehman representation of the spectral function^. 

First let us focus on the physically relevant case of 
U > and N = 2. To compare the numerical data with 
the weak coupling results, we have calculated the impu- 
rity density of states for different values of the interac- 
tion strength U. The results are shown in Fig. 2J The 
numerical data validates the weak coupling results for 
U/D < 0.3. In our NRG calculation we considered a flat 
band with constant density of states go = 1/2D, where 
D stands for the half bandwidth. In the lower panel 
of Fig. [4] the renormalized value of the hybridization, 
V ran is shown as a function of the interaction strength 
U. In NRG calculation, we have defined V lcn from the 
finite size spectrum directly. The finite size spectrum 
as a function of iteration number crosses over from the 
initial fixed point to the strong coupling one character- 
ized by single particle phase shifts S = tt/2 at around 
M*. M* is determined by the renormalized hybridiza- 
tion, A rcn = nV^go - A~ M */ 2 . We take M = M* 
when the energy of the first excited state exceeds 90% of 
its fixed point value. 

To answer the question whether an intermediate max- 
imum appears outside the weak coupling limit or not, we 
have performed calculation with very large values of the 
interaction strength up to U/D = 5.0. The results are 
shown in Fig. [S] 
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FIG. 5: (Color online) Comparison between weak coupling 
RG and NRG approach: The renormalized value of the hy- 
bridization, V ron as a function of the interaction strength U 
for N = 2 and N = 4. 



Our conclusion is that even for N = 2 such a non- 
monotonic behavior is found but the position of the max- 
imum as well as the shape of the curve for large U differs 
essentially from those obtained by weak coupling calcu- 
lation. It still remains a question whether for case of 
many channels the weak coupling calculation is reliable 
or not. To treat many channels with NRG is very chal- 
lenging but to see the tendency with increasing channel 
number, we performed the numerical analysis of the case 
N = 4. The results are plotted in Fig. O as well. Our 
data suggests that already for N = 4 the position of the 
turning point as well as the decay of the curve at large 
U is reproduced by the weak coupling calculation with a 
much better accuracy than in case of N = 2. 



IV. ANDERSON- YUVAL APPROACH 

In most of the physical cases the Coulomb interaction 
U dominates over the hopping term Vo- To overcome that 
difficulty, Anderson and Yuval^ introduced a path inte- 
gral method for the Kondo problem where the interaction 
U is described in terms of phase shifts while the hopping 
is treated as perturbation. The similarity between the 
Kondo and the present problem can be exploited in the 
following way: The complex time axis is divided into in- 
tervals and as it is shown in Fig. [5] where the solid lines 
represents the time interval when the impurity level is 
occupied and the light ones stand for unoccupied level. 
The conduction electrons can join the time line at the end 
points of the intervals where hopping Vq takes place while 
they can touch the time line at any other points due to 
the Coulomb interaction. Those are indicated by dashed 
lines which are labeled according to the channel indices 
oti . Thus the incoming and outgoing conduction electron 
lines should be connected all possible ways and finally a 
summation over all possible configuration of channel la- 
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FIG. 6: The solid lines represents the time interval when the 
impurity level is occupied and the light ones stand for un- 
occupied levels. The conduction electrons can join the time 
line at the end points of the intervals where hopping V a takes 
place while they can touch the time line at any points due to 
the Coulomb interaction. Those are indicated by dashed lines 
which are labeled according to the channel indices oti. 



bellings a, must be carried out. Some of the connections 
are indicated in Fig. [6] by dotted lines. The final result 
for the partition function can be given in analytical form 
as in Rcfs.^ili where the Kondo or the two level system 
problems were treated. The partition function has the 
form of a one dimensional Coulomb gas with appropri- 
ately defined vector charges C a (a = 1 . . .N — 1). The 
scaling equations are derived by eliminating short time 
intervals at the short time cutoff r with its initial value 
taken as the inverse bandwidth tq ~ D^ 1 . The inter- 
action Vq must be also made dimensionless by a factor 
~ r 1/2 

The phase shifts for electrons in case of filled and 
empty impurity levels are: 6 a — — arctan(go7rf ) an d 
S' a = — a.rcta.n(g TrU'), respectively. Only their differ- 
ence will appear in the scaling: z a — (S a — S' a )/ir. The 
phase shifts are limited: \5 a \, \S' a \ < tt/2. The Friedel 
sum rule requires J2 a z <* = ~ 1 which expresses that the 
one electron difference between the filled and empty sites 
must be screened by charge oscillations formed by the 
conduction electrons. The hybridization can be associ- 
ated with the fugacity 



V = V(g T a ) 1/2 cosSo 



(18) 



The interaction can be represented by charges at the in- 
teraction points as 



= ±{z a + 5 aQ ) 



(19) 



where the index ± labels hopping in and out of the impu- 
rity level assisted by an electron annihilation or creation 
in channel 0. 

The lengthy derivation of the scaling equations closely 
follows Refs. Ulil) and the final result is 



dy 
dlnr 

dC a 
dlnr 




-2y 2 C a 



(20) 
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or expressed in terms of phase z ( phase shifts) they are 



dy 
g?1iit 

dz a 
dliir 



= y 



[\-*>-\n% 

\ a / 



2{S a + z a )y 2 



(21) 



Here the scaling is carried out by increasing t to reduce 
the electronic bandwidth. The term ^y in the first line of 

Eq. (f2Tj) originates in the explicite factor r in the def- 
inition of y (|18[) and disappears from the corresponding 
scaling equation of V. 

The system of equations Eqs. (f2~Tj) must be solved for 
initial values y(r ) <gi 1 but z a can be arbitrary. The 
fugacity can either increase or decrease exponentially de- 
pending on the quantity (%+ \ E Q z a)- Similar expres- 
sions were obtained in Ref|7| by matching the perturba- 
tive results with expression in terms of phase shifts. 

The regions for attractive interaction and large enough 
repulsive one will be treated separately. In the first case 
(-zo - ^J2 a z l) < ( U < °> z o > 0). The solution is 
y/yo = ( r / T o)~ z °~^ ^= z " a thus y is decreasing, therefore 
z a {a = 0, . . . ) are slowly varying thus the r dependence 
can be ignored. Thus 



To 



(22) 



complete but at least it could be argued that it is inside 
the accuracy of the scaling equation. The weak coupling 
scaling result is very poor as it was expected (see FigH]) . 
For N — 4 all the methods give very similar results. 

The general solution of the scaling equations can be 
searched in the form C Q (r) = (C Q )i n itiaiC( r )- Then 



dlnr 
dy 



-2y 2 C 



1 



1 



dlnr \ 2 
The scaling trajectories are 



c 2 £c 2 



(23) 



% 2 (r)-C 2 (r)E^+41nC(r) =4^(0)-^^ . (24) 



During the renormalization y is fast increasing and the 
scaling is stopped where it reaches unity. Meanwhile £ 
decreases slowly and its renormalized value can be ex- 
tracted from Eq. ([M|) . The result in leading order in y 
reads as 



lnC(r) 



y 2 (r) 



y 2 



2 <t-->a 



(25) 



Outside that region the long time approximation for the 
conduction electron Green's function cannot be applied. 



The situation is different for repulsive interaction (z < 
0). There are two regions in that case. In the first one 
(— z — h E Q z 2 ) > and y is increasing. Then Eq. (|22]) 
is valid as far as y <C 1 is satisfied. There is, however, a 
crossover where (— zo~ \ ^ Q z 2 ) = to the second region 
where y decreases again and the screening dominates. 
The larger N the stronger the decrease is. 

The crossover between the increasing and decreasing 
regions is at 



As z a -s are very slightly renormalized thus the unrenor- 
malized values can be used and 5 a is independent of a 
(5 < 0). Thus the crossover is at z a — 2/N and then for 



N = 2 
N = 4 : 5 



5 = £ 

u 2 



QoU 
QoU 



oo 
1 



Comparing with the results of NRG in case N = 2 the 
turning point is at goU — ► oo thus the agreement is not 



V. WEAK COUPLING APPROACH FOR OUT 
OF EQUILIBRIUM 

Considering theoretical methods two ways can be fol- 
lowed: the Keldysh Green's function method or the 
calculation of the scattering amplitude by time-ordered 
perturbation theory. Here the second method will be 
followed, where the initial conduction electron states 
can be arbitrary non-equilibrium states and for the in- 
termediate and final states the actual non-equilibrium 
distributions are taken into account. That method 
has been earlier applied in the leading logarithmic 



approximation 



16.17.18 



, which is a generalization of An- 



derson's poor man scaling^. 

Here the extension of that method is presented to next 
leading order. For equilibrium first the Kondo model was 
treated by that way 2 ^. The basic idea is to calculate the 
development of the initial \i) state to the final |/) one, 
but in second order the renormalization of the norms of 
those states must be corrected also. Thus the scattering 
matrix clement to be considered is 



OO / 



n=0 



[{/I E rr^o^int \f)(i\ E rri^ifint 



(26) 
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(i) (ii) (iii) (iv) 




(v) (vi) (vii) 



FIG. 7: Panels (i)-(vi): The diagrams up to ~ U 3 order con- 
tributing to the numerator of Eq. (|26p . The diagrams should 
be decorated by the direction of the lines in all possible ways. 
The diagram of the self-energy is shown in panel (vii). 

where cu is the initial energy of state \i) and the Hamil- 
tonian is split as H = Hq + H ln t . In the present problem 
the correction of the normalization occurs for the impu- 
rity electron states while in the Kondo model for the spin 
states. Those normalizations appeared in the previous 
treatment as d 1 / 2 (cj) in Eq. (fT2|) . 

The scaling equation can be obtained after changing 
the cutoff D — > D — SD by adjusting the coupling con- 
stant to keep Tfi invariant for appropriate \i) and |/). 
In the following we use the original left and right states 
(q = L, R) with Vl = Vr as a start. In order to de- 
rive the scaling for coupling U the initial and final states 
should be atcffjO) where |0) is the non-equilibrium state 
at applied voltage eV, for which disregarding V the state 
is the noninteracting ground state. Considering the occu- 
pation rid of the d-level the occupation probability in the 
steady state must be determined in presence of eV. That 

I 



where the first correction is due to the vertex depicted 
in Fig[21 while the second one arises from the self-energy 
on the leg of diagram reduced by the denominator by a 
factor 1/2. This result agrees with Eq.lpB]) (N = 2) at 
eV = 0. Here taking the special case e& = the voltage 
eV serves as a low energy cutoff. 

As it has been mentioned earlier considering the d- 
level there is a steady state occupation rid- That value 
is determined from the ballance of the in and outflow of 
the conduction electrons. To determine it for Ed ^0 two 
other quantities must be known, namely the changes in 
the level position and the spectral function of the d-\eve\, 
id and Qd{s) due to the applied voltage. That calculation 
can be carried out numerically in a self-consistent way. 

The probability of scattering of an electron coming 
from the left (L) or right (R) into the d-\evel is denoted 



value will be (rid) — 1/2 for Ed = but in the general 
case Ed ^ it can depend on eV. The diagrams of the 
numerator up to ~ U 3 order are shown in Fig[7] (i)-(vi) 
where the diagrams should be decorated by the direction 
of the lines in all possible ways. The diagram of the self- 
energy is shown in Fig|7jvii). In logarithmic approxima- 
tion only the diagrams linear in log ^ are contributing to 
the scaling equations thus the relevant vertex corrections 
are (ii),(v) and (vi) while (iii) and (iv) are not as these 
provide log 2 The type (ii) diagrams with the parallel 
and antiparallel lines cancel each other. As the logarith- 
mic terms in diagrams (v),(vi) and (vii) come from closed 
electron loops which are independent of the actual value 
of fiL and fiR thus to the logarithmic term the left and 
right contacts contribute separately which should be in- 
dependent of the applied voltage. That simplification is 
not sustained in higher order contributions where the left 
and right lines simultaneously occur. The self-energy cor- 
rection in (vi) contributes by adding that to either of the 
incoming and outgoing d-lines. One of those corrections 
is cancelled by the denominator in Eq. (|2T)|) . As it is well 
known from the spinless fermionic case e.g. the X-ray ab- 
sorption proble m 11 ^ 3 the remaining diagram is cancelled 
by (v). Thus single logarithmic term does not remain. 
That is similar to the equilibrium case, see Eq.Q thus 
the invariant coupling J7i nv = const. (For the details see 
Appendix A.) 

In the following the renormalization of the hybridiza- 
tion depicted in FigJTJa and [2] is crucial where e.g. \i) — 
tft|0) and |/) = 4eJ°>- Keeping terms up to - VU 2 . Af- 
ter taking the denominator in Eq. (|2l)|) into account the 
final form of Tfi is 



(27) 

I 

by W£, R while the opposite process by W£i R . These 
quantities are 

W+ R = (1 - n d )2n go J Vl /R (E)Qd(E)f L/R (E)dE , (28) 
and 

W~ /R = n d 2TT Qo J Vl /R {e) Qd {e){l - f L / R (e))de , (29) 

where Vl/r(e) are determined from renormalization 
group equations with the appropriate infrared cutoffs and 
Il/r{ £ ) = f( £ i e ^/2) is the Fermi distribution function 
for the leads in presence of the voltage. Those will be 
taken at zero temperature T = 0. 



{Q\a ka Td)\0) = V a l + C/golog- 



D-E d 



eV 

a^f- - £ d 



1 9 9 D 
~U 2 Q 2 log- 
Z UJ 



Ed 
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The steady state is determined by This equation combined with Eqs. (f28l)([29|) gives 

jn d = W£ + W+ - Wl - Wx = . (30) 
I 



[ g d (e) [Vj(e)f(e + eV/2) + V%(e)f(e - eV/2)] de 

Jg d (e)[Vl(s) + Vl(e)]ds ' [ " } 



If electron-hole symmetry holds thus e d — 0, then n d — 
1/2. The next step is to determine the self-energy of the 
d-electron. The d-electron propagator is 

(oidtfx—L— # ld t| ) 

which can be simply developed because the occupied d- 
I 



level determines the time flow. The self-energy correc- 
tions appear also in the normalization. The effect of 
hybridization is just to give an extra broadening of the 
<i-level to be considered later. Without hybridization the 
self-energy is 



D+fi D+fi 

nu + i S) = U 2 gt J dt? J <%"(1 - /(O)Z(O u + i ,_^,_ Ed + i6 ( 32 ) 

-D+ii -D+n 



which can be evaluated as 



ReS(w) = U' a Qq I \oj - s d \ log ]UJ J dl + 2D log 2 ) , 



D 



and 



ImS(w) = ttU 2 q 2 {uj - e d )Q{uj - e d ) 



(33) 



(34) 



In the equilibrium calculation the term proportional to ~ 
w is contributing to the function d(ui) in Eqs.© and the 
last constant term ~ 2D log 2 is eliminated by the applied 
counter term to keep e d unrenormalizcd while ImS(w) 
is a Korringa type of relaxation. To be noted that the 
voltage does not occur as the energy goes directly into 
the electron-hole creation of the same electrode. As we 
mentioned at the end of Sec|TlJ that broadening is less 
important at small energies. The hybridization of the d- 
level gives the essential part of the broadening just like 



in the Anderson impurity model 

r(e)^2ng (Vl(e) + Vi(e)) 



(35) 



where the voltage dependent hybridization strength must 
be used. 

The the d-electron spectral function is 



1 



I»/2 



(36) 



7T (w - e d y + (r( w )/2)a • 

With the help of these quantities we are ready to calcu- 
late the current through th impurity: 

J=^(N R -N L ) = W^ + W+-W+-W£. (37) 

Combining this equation with the expression of W^, R 
given in Eqs. (f25))(|2l?|) the current takes the form 



n d 27r eo / [V£(e) - Vl(e)] g d {e)de - 2kq / [f(e - eV/2)V 2 (e) - f(e + eV/2)Vl(e)] g d (e)de 



The numerical calculation goes as follows: for a fixed value of eV we discretize the energy interval ui% G [— D + 
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U/D=0 

U/D=0.005 

U/D=0.03 

U/D=0.05 

U/D=0.10 



spectral function which scale as Qd{u = 0) ~ l/r(eV, U). 
For larger values of the voltage eV > A the d-level is 
experienced not only at the peak of the spectral func- 
tion (go(to = 0)) thus the current is not linear in eV any 
more. For eV 3> A the complete d- level contributes to 
the current thus it saturates and the asymptotic value is 
proportional to T(eV, U). 



VI. CONCLUSION 



|K , 1 , 1 , 1 , 1 

0.04 0.08 0.12 0.16 

eV/D 



FIG. 8: (Color online) Current obtained by weak coupling 
RG for A/D = 0.1, and U/D = 0,0.005,0.03,0.05,0.10. For 
this range of interaction strength the weak coupling method 
was reliable in equilibrium. 
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FIG. 9: (Color online) Current obtained by the scat- 
tering formalism for A/D = 0.1, and U/D = 
0, 0.10, 0.20, 0.30, 0.50, 0.60, 0.85. Note that the weak coupling 
result is reliable for U/D < 0.2 only. 



fi, D + fx] and calculate the renormalized hybridization 
VL/n(eV, uii) and the impurity self-energy. The latter 
is evaluated in such a way that the renormalized d-level 
position zero id = 0. By performing a sum over we cab 
calculate at the level occupation nd(eV) and the current 
I{eV) for the given value of the voltage. The result is 
shown in Fig[5]for U/D = 0, . . . , 0.10. 

In that regime of U the weak coupling RG gave good 
results in equilibrium (presented in Sec. lIII[) therefore one 
expects reliable results in out-of-equilibrium as well. As 
it is shown in FigJSl the current increases with increas- 
ing interaction strength. The results can be interpreted 
as follows: with increasing U the impurity spectral func- 
tion gets broadened since the hybridization is enhanced. 
In the linear regime (eV < A) the [/-dependence drops 
out since the increase of coupling to the leads r(eV, U) 
is cancelled by the decrease of the height of the d-level 



The resonant level model studied has very different 
behavior for attractive and repulsive interaction. This 
difference can be understood using the site representa- 
tion for conduction electrons in the strong U limit by the 
following argument: 

(i) in case of attractive interaction the particle on the 
d-level attracts electrons to pile up around the impurity 
occupying the next site thus the electron on the d-level 
is blocked for hopping to the conduction band. 

(ii) in case of repulsive interaction the site next to the 
occupied d-level is empty thus that electron can easily 
jump to the conduction band. 

All the methods predict that increasing the strength 
of an attractive interaction the hopping rate V is re- 
duced and for strong enough coupling it even goes to zero 
(see FigJ5]for Uqq < —0.25). The effect of orthogonality 
catastrophe reducing the hopping is less relevant because 
that have been already reduced by the vertex correction. 

In the repulsive case for large U the orthogonality 
catastrophe (self-energy correction) is reducing essen- 
tially the hopping rate. Thus the effective hopping V e g 
is first enhanced by the effect described above and can 
reach a maximum which is followed by a reduction due to 
the orthogonality catastrophe. The position of the max- 
imum can be pushed to lower energies by increasing the 
number of the screening channels N thus the perturba- 
tive result becomes more and more reliable. In case of 
N = 2 the latter method leads to a pronounced maxi- 
mum but the NRG indicates only slowly varying bump. 
In the Anderson- Yuval approach the maximum is even 
pushed to infinity. Thus the existence of the maximum 
is supported only by the NRG. 

Considering the time ordered scattering formalism the 
results are in accordance with the expectation of the weak 
coupling result for N = 2. The increasing U results in 
increasing current as first V is increased but for larger U 
due to the orthogonality catastrophe the current is essen- 
tially reduced (see Fig[S]). As the NRG does not give a 
sharp crossover for the hopping rate V r thus the corre- 
sponding effect in the current must be less pronounced in 
reality. Of course, for N ^> 2 the crossover must clearly 
exist. Unfortunately, in the work of Mehta and Andrei 
the crossover is in the range of reduced current, which 
is just the opposite what is expected on the ground of 
the physical picture established and results obtained by 
different methods for the hopping rate. 

Very detailed further studies of the Bethe Ansatz 
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FIG. 10: (Color online) The two relevant diagrams up to ~ (7 3 
contributing to the numerator of Eq. (|26[) . A Ward identity en- 
sures the cancellation of these diagrams for small frequencies 
meaning that no logarithmic term survives in second order 
thus the invariant coupling (7i nv = const. 

method are needed to understand and resolve the origin 
of the presented discrepancies. 
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APPENDIX A: CANCELLATION OF THE 

LOGARITHMIC TERMS IN THE 
RENORMALIZATION OF THE COULOMB 
INTERACTION 

The diagrams of the numerator of Eq. (|2"rj)) up to ~ U 3 
order are shown in Fig[7| (i)-(vi) and the diagram of 
the self-energy is shown in FigJTJvii). As noted earlier, 
the logarithmic terms come from diagrams (v),(vi) and 
(vii). Note that the self-energy correction in (vi) (See 
also FigfTUb) contributes by adding that to either of the 
incoming or outgoing d-lines. One of those corrections 
is cancelled by the diagram shown in Fig|7jvii)) in the 
denominator in Eq. (|26p . Therefore the two relevant dia- 
grams are those depicted in Fig llOl The contribution of 
FigfTUb can be written as 



D+fj, 

U 2 qI f # / — /(£"))/(£') — TTr~T7/ • ( A1 ) 

w - £d J J u + £' - 1," - £d 

-D+n -D+p, 



To get the purely logarithmic term we can now expand 
the integral in linear order in ~ (u> — Ed) and get the form 

I 



D+/J, D+fi 

uhi J # J d t"( l ~f(n)m') {Lj+e _\„_ £d)2 - (A2) 



On the other hand, the contribution of FigfTUb for 
small frequencies can be evaluated and the d-electron 
lines indicated by arrows occur twice just like the denom- 
inator squared in Eq. (|A2[) . That Ward identity ensures 



the cancellation of diagram shown in Fig[7Kv) by one of 
the leg one in Fig[7Jvi). This means that no logarithmic 
term survives in second order thus the invariant coupling 
U- lm = const. 
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